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Abstract. We study the asymptotic distribution of the number Zn of 
zeros of random trigonometric polynomials of degree N as N oo. It is 
known that as A'' grows to infinity, the expected number of the zeros is 
asymptotic to -^-N. The asymptotic form of the variance was predicted 
by Bogomolny, Bohigas and Leboeuf to be cN for some c > 0. We prove 
that ^w-EZjv converges to the standard Gaussian. In addition, we find 
that the analogous result is applicable for the number of zeros in short 
intervals. 



1. Introduction 

The distribution of zeros of random functions for various ensembles is 
one of the most studied problems. Of the most significant and important 
among those is the ensemble of random trigonometric polynomials, as the 
distribution of it zeros occurs in a wide range of problems in science and 
engineering, such as nuclear physics (in particular, random matrix theory), 
statistical mechanics, quantum mechanics, theory of noise etc. 

1.1. Background. Understanding the distribution of zeros of random func- 
tions was first pursued by Littlewood and Offord [LUlj . |L02j and |LU3j . 
They considered, in particular, the distribution of the number of real roots 
of polynomials 

(1) Pn{x) = ao + aix + . . . + a^x^ , 

of degree N with random coefficients an, as N — )■ oo. For the coefficients 
Un taking each of the values 1,-1 with equal probability 1/2, they showed 
that the number Zp^^ of zeros of Pn{x) satisfies 

(2) Zp^~-lniV 

vr 

for (l - O7v^oo(l))2^ of the vectors a G {±1}^. Later, Erdos and Of- 
ford |EOj refined their estimate. 

Kac [K] proved that the expected number of zeros Zp^ of the random 
polynomials ([1]) of degree A'^, this time a„ being Gaussian i.i.d. with mean 
and variance 1, is asymptotic to the same expression ([2]). His result was 
generalized by Ibragimov and Maslova |IMlj and |IM2j . who treated any 
distributions of the coefficients a„ , provided that they belong to the domain 
of attraction of the normal law: if each Ea„ = then the expectation is 
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again asymptotic to ([2]), though, if Ea„ 7^ one expects half as many zeros 
as in the previous case, that is 

EZp^ ~ ilniV. 
vr 

Maslova (Mlj also established the only heretofore known asymptotics for 
the variance of the number of real zeros Z, 

VarZp^ ~ -(1 - -) • IniV 
vr vr 

for the ensemble ([T]) of random functions. In her next paper [M2j . she went 
further to establish an even more striking result, proving the normal limiting 
distribution for Zpj^, as N ^ 00. 

The case of random trigonometric polynomials was considered by Dun- 
nage [DN]. Let Tat : [0, 2-k] ^ M be defined by 

N 

(3) T]\f{t) = a„ cos nt, 

n=l 

where a„ are standard Gaussian i.i.d, and Ztj^ be the number of zeros of 
T]\f on [0, 27r]. Dunnage proved that as — )• 00, E,Ztj^ is asymptotic to 

2 

and, moreover, that the distribution is concentrated around the expectation 
in the same sense as Littlewood and Offord mentioned earlier. 

The variance of the zeros for ([3]) was shown by Farahmand to be 0{N^/'^) 
in [F]. This estimate implies that the distribution of Ztj^ concentrates 
around the mean. 

Quails [QJ considered a slightly different class of trigonometric polynomi- 
als, 

N 

XN{t) = —1= / ( CLn sin nt + 6„ cos nt 



EZtj, 7^N, 



I f 

n=l ^ 



Let Zxpf be the number of the zeros of X]\f on [0, 2-k]. Applying the theory 
of stationary processes on X]\f, one finds that 



(iV + l)(2iV + l) 2 
^Zx^ = 2xl -^N, 



similar to ([3]). Quails proved that 



Zxj, - EZx 

N 



< C ■ iV3/4 



for some C > with probability 1 — OAr_>oo(l)- 

Bogomolny, Bohigas and Leboeuf [BEL] argued that the variance of Zxt, 
satisfies 

Var(Zx^) ~ciV, 
as N ^ 00, where c is a positive constant approximated by 
(4) c f« 0.55826 
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(this is equivalent to the formula (3.34) in |BBLj and the numeric value of 
A ~ 0.44733 immediately afterwards; one should bear in mind that they 
normalize the random variable to have unit expectancy). 

In this paper we study the distribution of the random variable Zxj^ in 
more detail. We will find the asymptotics of the variance Var(Zxjv) as well 
as prove the central limit theorem for the distribution of Zxj^ (see section 
II. 2p . We guess, but have not proved, that the same result may be true for 
Dunnage's ensemble ([3]). 

The zeros of random complex analytic functions were examined in a se- 
ries of papers by Sodin-Tsirelson (see e.g. [ST] ), and Shiffman-Zelditch (see 

e.g. m)- 

1.2. Statement of results. Let Xjv : [0, 27r] — )• M be Quails' ensemble 

1 ^ 

(5) Xjv(t) = a„ sinnt + fe„ cos nt, 



n=l 

where a„ and 6„ are standard Gaussian i.i.d. 

As usual, given a random variable Y, we define E(y) to be the expectation 
of Y. For example, for any fixed t £ [0, 27r] and A^, one has 

E(XN{tf) = 1. 

Above we noted Quails' result that 

(6) KiZx^) = 2^/x'2, 

where 

N 



n=l 



Ao := — > n 



N 6 

n=l 

We prove the central limit theorem for ^Xjv- 
Theorem 1.1. There exists a constant c > such that the distribution of 

Zx^ - EZx 

JV 

^/cN 

converges weakly to the standard Gaussian N{0, 1). The variance is asymp- 
totic to 

(7) Var(Zx^)~cAr, 

as N ^ oo, as predicted by Bogomolny-Bohigas-Leboeuf. 

We can compute the value of the constant c in Theorem 11.11 as 

4 2 

^ = ^^° + 7!' 

with 

oo 

/ON / [(1-5(2;)^)- 35' (X)2 /- —2 . 

(8) Co = / , ,^,0/0 (Vl-^* +-R arcsmi?)- 

j I {i-g{xyy/^ 


where we denote 



(9) 
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and 

g"{x){l-g{xf)+g{x)g'{xf 



More generally, for < a < 6 < 27r one defines Zxj^ (a, h) to be the number 
of zeros of Xn on the subinterval [a, 6] C [0, 27r]. It is easy to generalize the 
computation of the expectation ([6]) for this case as 



IT 

A priori, it seems that the behaviour of the number of zeros of in 
short intervals [aiVj^Af], shrinking as — )• cxd, should be more erratic than 
on the full interval. Surprisingly, just as in the previous case, we are able 
to find a precise asymptotics for the variance YaiZ^^aN ^bx)-, and prove a 
central limit theorem, provided that [aN-,bN] does not shrink too rapidly. 
We have the following Theorem: 

Theorem 1.2. Let < < < 2tt he any sequences of numbers with 
N ■ (6jv — cln) — ^ oo. Then as N ^ oo, 

Var(Zx^(aAf, 6Ar)) ~ c • "^-^ iV, 

zvr 

where c is the same constant as in Theorem \l.l[ Moreover, 

converges weakly to the standard Gaussian A^(0, 1). 

The proof of Theorem 11.21 is identical to the proof of Theorem II. H and in 
this paper we will give only the proof of Theorem 11.11 



1.3. Plan of the paper. To prove the central limit theorem we will first 
need to establish the asymptotic form ([7]) for the variance. This is done 
throughout sections [2] and [S] in section [2] we develop an integral formula for 
the second moment of the number of zeros, and in section [3] we exploit it to 
study the asymptotics of the variance. 

Sections H] and [5] are dedicated to the proof of the main statement of 
Theorem ll.il that is the central limit theorem. While the proof is contained 
in section SI a certain result, required by the proof, is proven throughout 
section \5\ 

1.4. On the proof of Theorem ll.lL As an initial step for the central limit 
theorem, we will have to find the asymptotics ([7]) for the variance. This is 
done throughout sections [2] and [3l 

While computing the asymptotic of the variance of Zx^^ , we determined 
that the covariance function rx of Xx has a scaling limit r^oix) = 
which proved useful for the purpose of computing the asympotics. Rather 
than scaling rx^ one might consider scaling Xx- 

We realize, that the above should mean, that the distribution of Zx^ 
is intimately related to the distribution of the number Zx of the zeros on 
(roughly) [0, A^] of a certain Gaussian stationary process Y[x)^ defined on 
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the real line M, with covariance function r = r^o (see section H?T]) . Intuitively, 
this should follow, for example, from the approach of |GSj . see e.g. Theorem 
9.2.2, page 450. Unfortunately, this approach seems to be difficult to make 
rigorous, due to the different scales of the processes involved. 

The latter problem of the distribution of the number of the zeros (and 
various other functionals) on growing intervals is a classical problem in the 
theory of stochastic processes. Malevich [ML| and subsequently Cuzick |CZ] 
prove the central limit theorem for Z]\[, provided that r lies in some (rather 
wide) class of functions, which include Too. Their result was generalized in 
a series of papers by Slud (see e.g. |SLj ). and the two-dimensional case was 
treated by Kratz and Leon |KLj . 

We modify the proof of Malevich-Cuzick to suit our case. There are 
several marked differences between our case and theirs. In their work, one 
has to deal with growing sums of identically distributed (but by no means 
independent) random variables (which will be referred to as a linear system); 
to prove the central limit theorem one applies a result due to Diananda [DNj . 
In our case, we deal with triangular systems (to be defined), applying a 
theorem of Berk [BR]. For more details about the proof, see section HTSl 

1.5. Acknowledgements. The second author wishes to thank Zeev Rud- 
nick for suggesting the problem as well as his help and support while con- 
ducting the research. In addition, he wishes to thank Mikhail Sodin, Ildar A. 
Ibragimov, Par Kurlberg and losif Polterovich for many fruitful discussions. 
We are grateful to Phil Sosoe for conducting some empirical experiments 
available in the longer version of this paper. We wish to thank Jonathan 
Keating for pointing out [BBLj . The authors wish to thank the anonymous 
referees for many useful comments and suggestions how to improve the pa- 
per. 

2. A FORMULA FOR THE SECOND MOMENT 
Proposition 2.1. We have 

(11) nzl,)-nZx,) = (\/r^+parcsinp^d^, 



where 

XX 1 2 {N + l){2N + l) 

(12) A2 = A2,7v = ^2^n = , 

n=l 



(13) 

and 
(14) 



N 



N ^ 

n=l 



COS nt 



2N 



sin {N + l/2)t - sint/2" 
sin (t/2) _ 



P = PNit) 



r"{l - 

A2(l 



r^) -I- (r')^r 



(r') 
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A similar but less explicit formula was obtained by Steinberg et al. [SSWZj . 
The rest of this section is dedicated to the proof of this result. 

The ensemble Xi\j is a centered stationary Gaussian process (meaning 
that the finite dimensional distributions are Gaussian with zero mean). An 
explicit computation with the double angle formula shows that its covariance 
function is 

rXr.itiM) = rxf,{t2 - ti), 
with the function on the right side as defined in 

Let / be an interval and X : 7 — )■ M be a mean zero stationary process with 
covariance function r. We assume that r(0) = 1 (i.e. X has unit variance) 
and furthermore that the sample functions of X are a.s. sufficiently smooth 
(e.g. twice differentiable) so that its sets of zeros is discrete. We have 

(15) \r{t)\ < 1 

for every t G / by the Cauchy- Schwartz inequality. We denote Z to be the 
number of zeros of X on I; Z is a random variable. 

In general, we have the following celebrated Kac-Rice formula (see e.g. |CL| 
for the expectation of Z: 

TT 

where |/| is the length of / (finite or infinite), and A2 = — r"(0). As men- 
tioned earlier, it was exploited by Quails to compute the expected number 
([U]) of zeros of trigonometric polynomials. 

In this section we find a formula for the second moment EZ^ of the 
number of zeros of any Gaussian stationary process X on I, assuming that 
its covariance function r is smooth. To determine E(Z^), we naturally 
encounter the distribution of the random vector 

(16) V = Vt,M ■■= {X{ti),X{t2),X'{h),X'{t2)). 

for some fixed ti,t2 G As an artifact of the stationarity of X, the dis- 
tribution of V depends only on t := t2 — ti. The covariance matrix of V 
is 



(17) S = S(t) 



/ 1 


r{t) 





r'{t) \ 


r{t) 


1 


-r\t) 








-r'{t) 


A2 


-r"{t) 


\r'{t) 





-r"{t) 


A2 / 



A B 
C 



The random vector V has a multivariate normal distribution with mean zero 
and covariance matrix E. For t G (0, 27r), S(t) is nonsingular (for X^ see [Q] 
and Remark 12.31). 



Lemma 2.2. Let X he a Gaussian stationary process, which almost surely 
has a continuous sample derivative such that the distribution of V is non- 
degenerate for ti ^ t2. Then 

E(Zi) - E{Zx) 

= jj dtidt2 jj \yi\\y2\(pt,,t2i0,0,yi,y2)dtidt2, 

[0,27r]x[0,27r] IR2 

where 4>ti,t2iui,U2,vi,V2) is the probability density function ofVt^^t2- 
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Remark 2.3. The original formulation of Lemma 12.21 from [CLj assumes 
that the Fourier transform of r has a continuous component, a stronger 
condition than stated. However, their proof works just as well in the less 
restrictive case as formulated above. Quails proved that the trigonometric 
polynomials ([5]) satisfy this assumption and thus we may apply Lemma 12.21 
to ([5]). Quails' argument can be generalized to higher moments: we use it 
to bound the third moment in Proposition lA.ll 

Remark 2.4. Let ipti,t2 be the probability density function of the random 
vector {X' (ti), X' {t2)) conditional on X{ti) = X{t2) = 0. Then we have 

(19) (pt,,t2{0,0,yi,y2) = - — ^===== 

27rY^l - r{t2 - hy 

(see also (j24p ). Therefore we may rewrite (jlSp as 

mz^^-mz )= ff ^[\x'ih)x'{t2)\\xih) = xfe) = o] dtidt2 

^ ^ J J V^-r{t2-hy 2n ■ 

Ixl 

We use this representation in the proof of Proposition 14. 3t as well as its 
analogue for the third moment in the proof of Proposition lA.ll 

We use Lemma 12.21 to derive the following. 

Corollary 2.5. Under the assumptions of Lemma \2JA one has 
(20) 

nZx^)-nZx^)= II _ ^2)3/2 (Vl-P' + Parcsmp)— 



Ixl 



where r = rx(t2 — ti), and p = px{i2 — ^i) with 

r"{t){l-r{tf)+r'{tfr{t) 
' A2(l -r(i)2) -r'(i)2 • 

Proof. Direct matrix multiplication confirms that 

i_/ (A-BC-^B^)-^ -A-^B{C - B^A-^B)- 
^ ~ \-C-^B\A- BC-^B*)-^ {C-B^A-^B)-^ 

so if n is the 2x2 "reduced covariance matrix", that is is the bottom 
right corner of then 

(21) n = C -B*A-^B. 

The matrix 0, is the covariance matrix of the random vector (X'(ti), X'(t2)) 
conditioned upon X{ti) = X(t2) = 0. 
Computing ()2ip explicitly, we have 

0, = /iJli, 



where 



with p given by and 



1 -p 
-p 1 



, , A2(l -r2) - (r')2 

(22) n := ^ % ^ ^ > 0. 
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Since is a covariance matrix, we have 

(23) \p\ < 1 

by the Cauchy-Schwartz inequahty. 
The easy to check identity 

yields 



I \o n 



(24) detS = det AdetO = (1 - r'^)iJ.'^{l - p^). 

Using (jlSp and the exphcit form of the Gaussian density (jji^ ,t2 ) we obtain 
(25) 

2n tf^7 ^ ff^,^, ff\ II | exp(-^j/l^~^j/*) dyidy2 

dtidt2 ff\ w \ f ^ -1 o-i t\ dyidy2 

mlml exp yil-^ y 



P 



Ml-r^)(l-p^)Vi -—--V 2- V (2vr)2 

M2 

U T r /"/"l II I 1 A dZldZ2 

zdtidt2 II |2i||z2|expl — z ' 



P 



V(l-r2)(l-/,2) ^^77'^"^' "V 2 ^ y (2^) 



where y = (yi, 2/2)) making the change of coordinates z = The inner 
integral is 



(26) 



|zi||z2|exp( - ^ + 2pziZ2 + ^i) ]dz 



4(1 — p^) ( 1 H ^ arcsin p 



computed by Bleher and Di [BDj . appendix A. Substituting this in ([25]) gives 
our result. 

□ 



We are finally in the position to prove Proposition 12.11 

Concluding the proof of Proposition \2.1\ We use Corollarv l2.5l on the trigono- 
metric polynomials X^- The integrand in ()20p depends only on t := t2 — 
(because of the stationarity of X^), which allows us to convert the double 
integral into a simple one. Proposition 12.11 then follows from the periodicity 
of the integrand. 

□ 



3. ASYMPTOTICS FOR THE VARIANCE 
Formulas (jlip and ^ imply the following formula for the variance 
(27) Var(Zx) = J + E(Zx), 



DISTRIBUTION OF ZEROS 



9 



where 



3.1. Evaluating the integral J in ([28]) . 
Proposition 3.1. 



dt 



(29) J=^Arfl + 0^^°^^ 



1/13 

37r " V ^ ' " V N~ 
where cq is defined by ([5]). 

Our key observation is that rxjv has a scaling Hmit, or, more precisely, 
we have 

/Ar(x) := rx„ ( — ) = h small error 

\m J X 

with m = N + ^ and x G [0,m7r] (treating [m7r, 2m7r] by symmetry), at 
least outside a small interval around the origin. It is therefore natural to 
change the integration variable in ()28p from mt to x, which will recover the 
asymptotics for J being linear with m, and thus also with N (see (j47p ). In 
fact, rather than introducing a new parameter m, it is also possible to use 
X = Nt; it results in a nastier computation. 

We will argue that it is possible, up to an admissible error, to replace 
the / = /at in the resulting integrand by g{x) := the latter being 

A'^- independent. The decay at infinity of the new integrand (i.e. with / 
replaced by g) imply that the integral will converge to a constant intimately 
related to ([8]). 

We divide the new domain of integration [0, vrm] (which is an artifact 
of changing the variable of integration x = mt; we also use the symmetry 
around t = it) into two ranges. Lemma 13.31 will bound the contribution of 
a small neighbourhood [0, 5] of the origin. The main contribution to the 
integral (|47p results from [6, vrm] , where Lemma 13.41 will allow us to replace 
/at in the integrand with A^-independent g{x) = 

Notation 3.2. In this manuscript we will use the notations A B and 
A = 0{B) interchangeably. 



Lemma 3.3. Let 



A'2(l-/(x)2)-/'(x) 



(30) M{x):= '\ ' ' {y/l-Rixr + R{x)aTcsmR{x)), 

(1 - f{xyy/^ 

where 

(31) f{^) = fN{x)=rxJ^) = (^-l^ ^ 



(32) R{x) = pN 
and 

1 + 

(33) A'2 = 



^ 2m -V 

2m 



x\_ f"{x){l-f{xf) + f{x)nx) 
m) A'2(l-/(x)2)-/'(x)2 
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There exists a universal 6q > 0, such that for any 

(34) 0<6 <6o, 

we have the following estimate 

5 

M{x)dx = 0{5^), 





where the constant involved in the "O" -notation is universal. 

Proof. First we note that 

(35) \Rix)\ < 1 

by the definition (|32|) of R{x) and (|23p . so that the definition of M(x) makes 
sense. 

We have to estimate f{x) and its derivative around the origin. Expanding 
/ and /' into Taylor polynomial around x = 0, we have 

fix) = l + amX^ + bmX^ + 0{x^), 

and 

fix) = 2amX + 4bmX^ + 0{x^). 

with 

_ (2m + l)(12m^-7) 

2880m3 -^^^>' 
and the constants in the 'O '-notation being universal. 
Thus, 

1 - /(x)2 = -2a^x2 - {2bm + al)x^ + ©(x^) > x^, 

and 

A'2(l - f{xf) - fixf = i±2k(-2a„x2 - {2bm + al)x^ + 0{x^)) 



8640m4 
Now 



4x^(a^ + Aambmx'^ + Oix"^)) 

'^^ + 0{x^) < x^ 



64m*^ + 24m^ - lOBm^ - 94m - 21 4 g 4 



1 ^ \/l — + y arcsin y < tt/2 



for every y G [—1,1], so combining the last three displayed equations, we 
obtain 

4 

M(x) < i^i- = X, 

X'^ 

and the Lemma follows. 

□ 
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Lemma 3.4. Let 

(36) 
and 



5 > (m/2)-i/9 
6 < X < nm. 



Denote 

Mi(x) :=M(x)-A'2, 
where M{x) is given by ()30p . Then, form sufficiently large, 

{l-gixf)-3g'ix)\ 



Mi{x) = - 



(1 - 5(x)2)3/2 



1 - + R*{x) arcsmii*(x)) - 1 



where, as usual, g[x) and R* are given by Q and (jlOp respectively. 

Proof. We will approximate / and its first couple of derivatives by g and its 
first couple of derivatives. 
For 5 < X < TTm, we have 



fix) 



(37) 



smx 



1 



^(1 + 0(S)) /2m-l 



1 



2m • sin(x) 



l + O 



1 



1 



2m - 1 



g{x) + O 



where to justify the second equality, we use the explicit coefficient of the 
second summand in the Taylor formula for the sine. By a straightforward 
computation of Taylor's formula, one easily obtains the sharp estimate 

„2 



fix) = g'{x) + O 



9'{x) 



1 x' 
— + 



+ 



g{x) 



m m- 

for < a; < vrm. Noting that g{x) = O ( p^^pri ) , d'ix) = O 



m^J 



we have 



(38) 



f'{x)=g'ix) + o(^)+o(— 
\m'^ J \mx 



for Q < X < vrm. For < x < 1 we have the better estimate 

f'{x)=g'{x) + 0(-). 

Similar statements hold for /"; in particular, for 5 < x < Tim we have 
(39) f"{x)=g"{x) + o(-^ + ^\ 



Next, we apply (|37j) to obtain 



l-/(x)2 = (l-5(x)2). 1 + 



52 



mx 



where we used for x > 



1 - g{xf > 1 - g(bf = ^ ^2 
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Thus 

(40) (1 - fixf) = (1 - gi.fr^l- • (l + O (^) ) : 

where we used the assumption (j36p to bound away from 1. 
By the above, we have 



(41) 



A'2(l - f{^f) - = A'2(l - 9{x?) - g'{xf + O {-^ 



since for x > 5, 
(42) 

A'2(l - g{xf ) - g'{xf > ^(1 - g{xf) - g'{xf > ^(1 - g{5f) - g'{5f » 5\ 
Next, we have 

f"{x){l - f{xf) + f{x)nxf = g"{x){l - g{xf) + g{x)g\xf + o{^ + 



mx 



using (p6|) again. Therefore, 

j.r .^ g"{x){l-g{xf)+g{x)g'{xf ( 1 x 
^ ' g{xf)-g'{xY \6^mx 6^m?)' 

exploiting ()36p once more as well as ()42p . 
Note that 

g"{x){l - g{xf ) + g{x)g'{xf g"{x){l - g{xf) + g{x)g'{xf 



- g{xY) - g'ixy 1(1 - 5(x)2) - g'ixy + 0(i) 

g"ix)il-gixr)+gix)g'ixr ^^f 1 



^{l — g{x)'^) — g'{x)'^ \6^mxJ' 

where we use twice as well as again. 
All in all we obtain 

(43) Rix) = R* (x) + O f + ^ 



^mx 

where R* is defined by (jlOp . It is important to notice that (j35p implies 

|i?*(x)| < 1, 

since for any fixed x, 

R{x) R*{x), 

as m — )• oo by and . 
Notice that 

(44) «-W = 0(sfe)- 
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again by dMI). Therefore for j; > 1/5"^, R*{x) is bounded away from 1 so 
that for a small e > one has 



y^l - {R* + e)^ = y^l - R*{x) ■ J 1 



2ei?*(x) + e2 



1 - 



= - R*{x) ■ - 0{eR*{x) + e2) = ^1 - + 0{eR*{x) + e^). 

It yields 

and similarly 

arcsini?(x) = R*{x) aicsin R* (x) + O 
Therefore for x ^ 1/5'*, we have 



- + arcsin i?( 



a; 



(45) /-— — ^ ... . 1 



1 + arcsin i?* (x) + ^+ „ , . 



On the other hand, for 5 < x ^ 
— R{xy + R{x) arcsin 



l-R*{xf + R*{x) arcsin ^* (^) + O ( ^ + ^ 

9 /II 

1 - i?*(2;)^ + i?*(x) arcsin + ^ttj h 75-^ 



since the derivative of the function 



X I—)- \/l — + X arcsin x, 

namely arcsin x, is bounded everywhere in [—1, 1]. Therefore (|45p is valid 
for X > 5. 

Also we have 



(46) A/l-^*(a;)Vi?*(x)arcsini?*(x) = l + 0(i?*(x)^) = l + O 



1 



5^x2 / ■ 



by m- 

Collecting and (00]), we have 
A-2(l-/(x)^)-/W ^ \'^{l-g{xf)-g'{xf ( ( 1 

(l-/(x)2)3/2 (l-5(x)2)3/2 \^ ^ \^<56^a; 

_l/3(l-5(x)2)-5'(x)2 J_ ^ 1 



(1 — (7(x)2)3/2 6771 y^Sj^^, 
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Together with (j45p it gives 

\2\ fit^\2 



1/3(1 -ff(x)^)-g^(:r) 

(1 - g(2;)2)3/2 



1 - i?*(x)^ + R*{x) axcsm.R*{x] 



{l-gixf)-3g'{xf 
(l-5(x)2)3/2 



1 - R*{x)'^ + R*{x) arcsin - 1 



by (06]). 



□ 



Proof of Proposition \3.1\ Noting that the integrand in (|28p is symmetric 
around t = tt, denoting m := N + ^ and changing the variable of integration 
mt to X in (|28p . we find that J is 
(47) 



nm 

J= — j y (1 _ /(^)2)3/2 \^/l-R{x?+R{x)&TCSUiR{x)y\^ 



dx. 



where /, R and A2 are defined in ([3T]) . ([32]) and ([33]) . 

We divide the interval into two ranges: Ii := [0, 6] and /2 = [S, Trm], for 
some parameter 6 = S{m) > 0. On Ii we employ Lemma [3.3l to bound (from 
above) the total contribution of the integrand, whereas we invoke Lemma 
13.41 to asymptotically estimate the integral on l2- The constant 6 has to 
satisfy the constraint of Lemma 13.41 namely ([36p . The constraint of Lemma 
3.31 ([M]) . is satisfied for m sufficiently large, provided that 6 vanishes with 
m. To bound the contribution of A2 to the integral on Ii, we use the trivial 
estimate A2 = 0(1). 

Hence we obtain 

(48) 

Trm, 

{l-g{xr)-3g'{xy 



Trm 

4m 
3^ 





/ 



(l-5(x)2)3/2 



l-R*{xf + R*{x) arcsin E*(j;) - 1 



^,logm 1 , , 
Note that for a large x we have 



(l-<7(x)2)-3ff'(x)2 = l + o(^) 
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and the definition (llOp implies 



so tliat 



I- R*{xf + R*ix)sxcs\nR*{x) = l + 0{R*ix)) = 1 + 0(4- 

Plugging the estimates above into the integrand of ()48p shows that the 
integrand is O(^). Hence 



{l-g{xf)-?>g'{xf 
(l-(7(x)2)3/2 



1 - R*{xf + R*{x)B.Tcs{nR*{x) - 1 



oo 

_ 4^ 
37r 



logiV 1 
+ 0{^ + ^) + 0{5N). 

We finally obtain the statement ([29|) of the present proposition upon choos- 
ing 

'1 Ar\ 1/13 
, log A* \ ' 

:= 



N 

□ 

3.2. Concluding the proof of the variance part ([7]) of Theorem II. IL 

Proof. Proposition 12.11 together with Proposition 13.11 and ()27p imply 

Var(Zx) = J + nZx) ~ — iV + ^iV = cN. 

Svr V3 

It then remains to show that c > 0. As mentioned earlier, Bogomolny- 
Bohigas-Lebeouf [BBL] estimated c ~ 0.55826 and it is possible to use nu- 
merical methods to rigorously obtain c > (see the longer version of the 
present paper). 

There exists a more systematic approach though. One can construct a 
Gaussian process y^(j;) on ffi with the covariance function Vaoi^^ 
(see section HTT|) . In this case, we denote again 



sin X 

X 



1 

A2,oo = -?'^(0) = -. 

For r > 0, let Zoo{T) be the number of the zeros of Yoo on [0,T]. 

By the general theory of stochastic processes developed in [CLj , one has 

T 

^Zoo{T) = — A2,oo- 
TT 

Using the same method we used to compute the variance of X^v, it is not 
difficult to see that 

VarZoo (T) ~ cT, 

where c is the same constant as in Theorem ll.il provided that c > 0. It was 
proved by Slud [SLJ , that it is indeed the case for a wide class of covariance 
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functions r(x) which contains our case r = Too- Moreover, Slud (following 
Malevich [MLj and Cuzick |CZj ). established the central limit theorem for 

Zoo(r)-EZoo(T) 

□ 

4. Proof of Theorem 11.11 

In this section we pursue the proof of the central limit theorem. The 
main probabilistic tool we use is a result of Berk [BRj . which establishes a 
central limit theorem for a triangular system of random variables defined 
below. We start however with a general remark about the zeros. 

4.1. A general remark about the distribution of the zeros. Let 

Yoo : M — 7- M be a Gaussian stationary process with the covariance function 
''00(2;) = Such a process exists, since Too has a nonnegative Fourier 

transform on R (this is related to Bochner's theorem, see e.g. |CL] . page 
126). Moreover, we may assume with no loss of generality that Y^o is al- 
most surely everywhere differentiable. In fact, one may construct Yqo using 
its spectral representation. Alternatively, one may use the Paley- Wiener 
construction Y{x) = ^ ■ 

To define all the processes on the same space, we will assume that Yj\[ are 
defined on R by periodicity. We have the convergence ry^(x) — >• ry^(x), 
as iV — 7- 00. This implies the convergence of all the finite-dimensional dis- 
tributions of Y/v to the finite-dimensional distributions of 1^. By theorem 
9.2.2 [GSJ, for any continuous functional (j) : C{[a,b]) — )■ M, the distribu- 
tion of (j){YN) converges to the distribution of (j){Yoo) (one may easily check 
the additional Lipshitz-like condition required by that theorem). Thus one 
could model a "generic" statistic of Atv by the corresponding statistic of Y^o 
on intervals, growing linearly with N. 

The convergence r^v — >• ?'oo suggests that the distribution of the number 
of zeros of A^r on the fixed interval [0, 27r] is intimately related to the distri- 
bution of the number of zeros of the fixed process on growing intervals. 
The particular case of the latter problem when the process is Gaussian sta- 
tionary (which is the case in this paper), has been studied extensively over 
the past decades. 

4.2. Triangular systems of random variables. Let C = {Ik} be a (finite 
or infinite) sequence of random variables (which we will refer as a linear 
system), K its length [K = 00 if £ is infinite), and M > 1 an integer. We 
say that C is M-dependent, if for every i, j > 1 with i — j > M, {1^ : k < i} 
and {/fe : k > j} are independent. For example, a 0-dependent linear system 
C is independent. 

One is usually interested in the distribution of the sums 

N 

Sn = ^ h 

k=l 
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as — )• oo, that is if -ftT = oo. In this situation one employs a special version 
of the CLT due to Diananda |DNj . For a process X{t) on M we denote Zx{T) 
to be the number of zeros of X on [0,T]. Cuzick |CZ] employed Diananda's 
result to prove the central limit theorem for Zx{T) as T — )• oo for a wide 
class of stationary processes X. A more basic version of this theorem was 
applied earlier by Malevich [ML) to obtain a similar, but more restrictive 
result. 

Diananda's result applies to finite sums of random variables of linear 
systems. The situation in our hands is somewhat different, namely, of a so- 
called triangular system (or array) of random variables. A triangular system 
of random variables is the correspondence K{N) : N — )• N U {oo}, together 
with a linear system 

Cn = {zN,k : l<k< K{N)} 

for each G N. We will use the notation Z = {z^^k} to denote a triangular 
system. 

Let M = M{N) be sequence of integers. We say that Z is M-dependent, 
if Cn is M(A^)-dependent for every A^. 

Given a triangular system Z, we are usually interested in the asymptotic 
distribution of the sums 

K(N) 

Sn = ^ ^AT.fc, 
fe=l 

as — 7- oo. Note that here, unlike the linear systems, both the number of 
the summands of Sm and the summands themselves depend on A^. We have 
the following theorem due to Berk [BR] , which establishes the asymptotic 
normality of Sn for M-dependent triangular systems: 

Theorem 4.1 (Berk [BR]). Let Z = {zN,k ■ I < k < K{N)} be a M- 
dependent triangular system of mean zero random variables. Assume fur- 
thermore, that 

(1) For some 5 > 0, Elzjy.fcP^*^ ^ ^i; where Ai > is a universal 
constant. 

(2) For every N and 1 < i < j < K{N), one has 

Vai{zN,i+i + . . . + ZNj) < (j - i)A2 

for some universal constant ^2 > 0. 

(3) The limit 

Var(zAr,i + ... + zn,k) 



lim 



Af->oo K 

exists and is nonzero. Denote the limit f > 0. 
(4) We have 

M = o(^K2^ 

Then 

zn,i + . . . + Zn,K 

converges weakly to A'^(0, 1). 
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Note that condition S] requires, in particular, that as N ^ oo, one has 
X — )• oo. Berk's result was recently generalized by Romano and Wolf |RWj . 

4.3. Plan of the proof of Theorem 11.11 We define the scaled processes 

Yn{x) := Xn 

on [0,27rm], where we reuse the notation m = N + 1/2 from the proof of 
Proposition 13.11 Let us denote their covariance function 

r^ix) = rvj^ix) = rXffi^^ = fN{x), 

with Jn defined by (f3T]) . It is obvious that 

the number of zeros of Yn{x) on 

I'n := [0,2™]. 

It will be sometimes more convenient for us to work with Yn{x) and rAr(x) 
being defined (by periodicity) on 

In := [— vrm, Trm] 

rather than on 

We are interested in the number Zyj^ of zeros of Yn on intervals In, whose 
length grows linearly with A^. Divide In into (roughly) subintervals lN,k 
of equal length with disjoint interiors (so that the probability of having a 
zero on an overlap is 0), and represent Zyjsr, almost surely, as a sum of 
random variables ZN,k, the number of zeros of Yn on lN,k- The stationarity 
of Yn implies that for a fixed N, ZN,k are identically distributed (but by no 
means independent). 

We, therefore, obtain a triangular system 

Z = {ZN,k = ZN,k — EZtv.A;} 

of mean zero random variables with growing rows. Just as ZN^ki the ran- 
dom variables ZN^k are identically distributed. We will easily show that Z 
satisfies the conditions [2]l3] of Theorem 14.11 Moreover, we will see later that 
condition [1] holds with 5=1 (see Prop osit ion IA.lt here we deal with a more 
complicated case of mollified random variables defined below; an easier ver- 
sion of the same argument applies in case of Z, however it does not give the 
CLT due to the lack of independence). 

The main obstacle to this approach is that the random variables ZN,k 
are not independent (and thus neither are ZN,k)- In fact, we may give an 
explicit expression for 

Cov{ZN,ki, Znm) 

in terms of an integral, which involves rjv and its derivatives. The station- 
arity of Yn implies that Cov{ZN,ki, ^N^k-i) depends on the difference k2 — ki 
only (that is, the discrete process ZN^k with A^ fixed is stationary, provided 
that the continuous process Y/v is). 

To overcome this obstacle, we notice that rjv(x) and its couple of deriva- 
tives are small outside a short interval around the origin, and, moreover. 
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their mass is concentrated around the origin. This means that the de- 
pendencies between the values and the derivatives of Yi^{x) on I^^kx and 
those on lN,k2 "small" for \ki — A;2 1 sufficiently large. Thus the system Z 
is "almost M-independent" , provided that M = M{N) is large enough (it 
is sufficient to take any sequence M{N) growing to infinity; see Proposition 

1131). 

One may then hope to exchange the process Y^r (and thus the system Z) 
with a process Y^^ (resp. Z^^), where M = M{N) is a growing parameter, 
so that the distributions of the number of zeros Z^v = Zyj^ and Z^^ = ZyM 

of Yn and Y^^ respectively, are asymptotically equivalent, and the above 
properties of the original system stay unimpaired. In addition, we require 
Z^ to be M-dependent (or, rather, const ■ M-dependent). To prove the 
asymptotic equivalence of Z]\f and Z^ we will evaluate the variance of the 
difference Yai{ZN — Z^) (see Proposition 14. 3|) . 

To define Y^^ , we introduce a function = ■ Sm, where \Sm\ < 
1 is a sufficiently smooth function supported on [—const ■ M, const ■ M] 
approximating the unity near the origin with a positive Fourier transform 
on the circle In (with end-points identified). We require that (and a 
couple of its derivatives) preserve 100% of the mass of r^v (resp. a couple 
of its derivatives) (see Lemma |5.1|) . We then construct Y^ with covariance 
function using some Fourier analysis on In- It is important to observe 
that the covariance function being supported essentially at [—M, M] means 
that Z'^^ is (roughly) M-independent, in the periodic sense. 

To get rid of the long-range dependencies, which occur as a result of 
the periodicity of Yn, we prove the central limit theorem for positive and 
negative zeros separately (see the proof of Proposition 14. 4|) . Namely we 
define Z*^'"*" (resp. Z*^'~) to be the number of zeros z of Y^^ with z > 
(resp. z < 0), and 

_ ^M,+ _|_ ^M,- 

almost surely. We are going to prove the asymptotic normality of the distri- 
bution of Z^^'"*" and similarly, of Z^^'^. We will prove that this will imply 
the asymptotic normality of the sum Z*^. 

Concerning the choice of M = M(N), on one hand, to well approximate 
Y/v, M has to grow to infinity with N. On the other hand, condition [5] 
of theorem 14.11 constrains the growth rate of M from above. The above 
considerations leave us a handy margin for M. 



4.4. Some conventions. In this section we will use some Fourier analysis 
with functions defined on the circle In := [— vrm, 7rm] (or equivalently, := 
[0,27rm]). We will adapt the following conventions. Let / : /^v — )• M be a 
real-valued function. For n G Z, we define 

fin) = [ me-'^^^. 
J \/2TTm. 
In 
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If / is a real valued even nice function, then 

r{x) = f(0) • ^= + y V2r{n) ■ """^IZ^ . 

^ n=l ^ 

With the above conventions, if f, g : In ^ M are two functions, then 

(/•5)(n) = -^(/*5)(n), 
\'2TTm 

and 

(/ * 9){n) = \/27rmf{n) ■ g{n). 
For the real valued even functions, the Parseval identity is 

oo 

11/11 W.) = ii/ii'^(z) = /W + E 2/»'- 

n=l 

4.5. Proof of Theorem II. IL 

Proof of Theorem \l.l[ Let Yn{x) = Xn{^), and for notational convenience, 
we assume by periodicity, that Yat and its covariance function are defined 
on Jtv := [— vrm, Trm]. One may rewrite the definition of Yat using 

1 °° 1 / \ 

(49) r7v(x) = fAr(O) • , + V \/2rAr(n) • —= cos ( — x ) , 

V27rm ^J'nm \m ) 



n=l 

as 

lV(x) = v^^^^^ao 

(50) 



+ 2^/'^v^fAr(n) • - — a^cos \ —x \ +6„sin \ 

^1 (7rm)V4^ \ra ) \m J J 

where a„ and 5„ are (0, 1) Gaussian i.i.d. One may compute tat to be 



otherwise 

It is easy to identify (|50p as the spectral form of Y/v on the circle, analo- 
gous to the well-known spectral theory on the real line (see e.g. |CL) . section 
7.5). The spectral representation proved itself as extremely useful and pow- 
erful while studying various properties of stationary processes. 

Let < M < vrm be a large parameter and X[-m,m] characteristic 
function of [-M,M] C Ij^. Define 

Sm{x) = , 

where (•)*' stands for convolving a function I times to itself, and the universal 
constant C > is chosen so that 5Af(0) = 1. The function Sm : /at — )• M is 
a piecewise polynomial of degree 7 in independent of M. It is a 6-times 
continuously differentiable function supported at [— 8M, 8M]. For |x| < 2M, 
for example, 

(51) 5„(.) = 1 + (|j)' + (^) ' + (s)" + '^(^, ' 
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for some constants 5i, . . . , ^4 G M, which may be easily computed. 
We define the moUified covariance function r^^ = : I]\f 

(52) r^{x):=rN{x)-SM{x), 



with the Fourier series given by 

1-1 

(53) rfj{n) = • (f^ * SM)in) = :rT-Axi<\n\<N * 5'm)(?^) > 0, 

yzirm ^^V 

since 

(54) SM{n) = ^^^^ • {X[-MM]i^)f > 0- 



One may compute exphcitly the Fourier transform of X[-m,m] to be 
(55) X[-M,M]("-) = 



n=0 

-■^ — sm , n 7^ 

The nonnegativity of Sm allows us to construct a process {x) on /^r 
with covariance function Tvm = r¥ as 

^^^^ 2V4 / (n \ ^ . (n 

+ L V • (^^ (;;^^J + 

the RHS being almost surely an absolutely convergent series, uniformly w.r.t. 

X. 



Remark 4.2. We assume that the and 6„ for n < in (j56p are the 
same as in ([50]) . so that Y^^(x) converges in to lAr(a;) (see Lemma [53]) . 

Let M = M{N) be any sequence of numbers growing to infinity, satisfying 
M = o{N'-'^). Proposition [Ql then implies that as 00, ^^^^ is 

asymptotically normal. Proposition 14.31 states that 

Var(Zi^^-Z,v) =o(Var(Z^)), 

so that the distribution of '^^^^^ ^ is asymptotically equivalent to that of 
; which implies the statement of Theorem II. li 

n 

Proposition 4.3. Suppose that as N ^ 00, we have M — )• 00. Then one 
has 

\aT{ZN - Zfl) = o{N). 
The proof of Proposition \4-S\ is given in section O 
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4.6. Proof of CLT for ZfJ . The main result of the present section is 
Proposition 14.41 which estabhshes the central limit theorem for the mollified 
random variable . 



Proposition 4.4. Suppose that as N ^ oo, M = o{N^ ). Then for 
N — )• oo, the random variables ^ ^ , weakly converge to the standard 
Gaussian N{0, 1). 

Proof of Proposition \4-4\ Recall that is the number of zeros of on 
Ii\f = [— 7rm, vrm]. We are going to prove a central limit theorem for the 
number of zeros on := [0,7rm] (denote similarly := [— vrm, 0]) only. 
We thereby denote zff'~^ (resp. Z^~) the number of zeros of Y^^ on 
(resp. /j^), and analogously, and Z^ will denote the number of zeros 
of Yat on I'^ and respectively. This also implies a central limit theorem 
for Z^'~ by the stationarity. The problem is that Z^'^ and Z^~ are not 
independent, so that writing Z^ = Z^'~^ + Z^'~ a.a. does not imply the 
asymptotic normality for the sum. Therefore we will have to come up with 
a more gentle argument in the end of this proof. 

Let -L > be any integer, which we will keep fixed throughout the proof. 
We divide I'^ into subintervals 



for 1 <k < LN, and denote Zfff^ the number of zeros ofY^'ix) on lN,k- 
Recall that, as a function on [—vrm, vrm], rfj is supported on 

[-8-M,8-M]. 

Therefore, if i — j > 8LM, the random variables {Z^j^ : k < i} are inde- 
pendent of {Z^j^ : k > j}. 

We apply Theorem l4.1l on the M = const • M-dependent triangular system 

Z'^ = {Zl,\, := Zl,\, - EZ*^, ■.l<k< K{N)} 

with K{N) = NL. Since with probability 1 neither of Yjf have zeros on 
the overlaps of In.u-, we have 

NL 
k=l 

almost surely, so that to finish the proof of this Proposition, it remains to 
check that Z^^ satisfies conditions [T][3] of Theorem 14.11 

Proposition lA.ll implies that condition [1] is satisfied with 5 = 1, provided 
that we choose L large enough. Since M ~ const -M and K{N) ~ const- 
the assumption M = o(N^/'^) of the present Proposition is equivalent to 
condition [H 

Condition El of Theorem 14.11 is equivalent to Yav{Z^'^) ~ ciN for some 
ci > 0. An application of ([7]) together with Proposition 14.31 and the triangle 
inequality, imply that 

Var(Zj^) ~ cN. 
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One may also derive the corresponding result for Z^'~^ , starting from Var(Z^) 
|iV (the proof follows along the same lines as the proof of ([7])) and using 
(|61|) with the triangle inequality. 

It then remains to check that satisfies condition [2] of Theorem 14.11 
Using the same approach we used in the course of the proof of ([7]) , one may 
find out that 
(57) 

Var(Z*^,+i + . 



+ Z 



M ) 
N,j) 



7rm 



r(x)^) 



r'{xY 



r\x) ) 



+ p[x) arcsin /9(x)) — A 



M ' 
2,N 



dx + ( j 



where we use the shortcuts r 



N 



(0), and 



p{x) = pfJix) 



r"{x){l — r(x)^) + r'(x)^r(a;) 
A'2(l -r(x)2) -r'(x)2 ' 



We have 



since ^ < 2 and A 



M ' 

2,N 



(j-^);^yA|57«(j-.), 

= 0(1). It remains therefore to bound the integral 



in (j57p . which we denote J. We write, denoting r := (j 



'') LN- 



(58) J«(j-i) 



j2 



.2^3/2 



■ ( yl — fP' + p arcsin p) — A 



M' 



dx. 



It will suffice then to prove that the latter integral is uniformly bounded. 
Let K^{x) be the integrand. Expanding Kf} [x) into Taylor polynomial 
around the origin, as we did in the course of the proof of ([7]) (see the proof 
of Lemma 13. 3p . we find that K^{x) is uniformly bounded on some fixed 
neighbourhood of the origin (say, on [0,c]). We claim, that outside [0,c], 
the integrand is rapidly decaying, uniformly with N. 

It is easy to see that r{x) being supported at [0, const ■ M] implies that 
K{x) is supported in [0, const ■ M] as well (note that we exploit here the 
fact that by counting only the positive zeros we disregard the dependencies 
between zeros near —nm and vrm). Moreover, on [c, const • M], {x)\ <^ 
where the constant involved in the "<C" -notation is universal. Therefore 
the integral on the RHS of (f58|) is uniformly bounded, so that J <^ (j — i), 
which verifies condition [2] of Berk's theorem. 

This concludes the proof of the asymptotic normality for Z^'~^ (and also 
Z^'~). Having that result in our hands, we define the random variables 
Z^'^ and Z^'~ to be the number of zeros of Y^^ on [8M, vrm - 8M] and 
[—vrm + 8M, —8M] respectively. The random variables Zj^ ' are indepen- 



dent, since rff is supported on [— 8M, 8M]. 



In addition, let zff'^ , Z^'^, Z^'g 



and Z 



Y^' on [0,8M], [vrm-8M,vrm], 



N,L 

8M, 0] and [-vrm 



be the number of zeros of 
vrm+8M] respectively. 
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We have 

VarZ^^^, \aiZ^'^, VarZ*^^, YcirZ^^ ^ ^ 
(59) ' iif ' 11 f ' 

= o(VarZ^^'+), o(VarZ^^''), o(VarZi;f) 

by condition of Theorem 14.11 (which we vahdated) . Therefore 



and 



yM- _ yM- _ yM,- _ yM- 



converge to the Gaussian distribution 



The independence of Z^'^ then implies the asymptotic normahty of 



^n'~^ + Z^' , and finahy we obtain the asymptotic normahty of 
again by ([MI)- 



n 



5. Proof of Proposition 14.31 



5.1. Introduction and the basic setting. Recall that we have the pro- 
cesses Y]\[{x) and Y^{x), defined on 1^ = [— vrm, vrm] and are interested in 
the distribution of Z^r and Z^^ , the number of zeros of Yat and Y^^ on 1^ 
respectively. The goal of the present section is to prove the bound 

(60) Var(Z^ - Zf^) = o{N) 



on the variance of the difference. For notational convenience, we will con 

(I'^Sp. Z^ 



sider only the positive zeros, that is, let Zir (resp. Z^'^) be the number of 



zeros of Yn (resp. Y^ ) on = [0, vrm]. We will prove that 

(61) Var(Z+-<'+)=o(iV), 
and by the stationarity, it will also imply 

(62) Var(Z^-Zj;f'") =o(iV), 

where we denoted the number of negative zeros in an analogous manner. 
Finally, (|6ip together with (|62p . will imply (j60p . by the triangle inequality. 

Let S > and i? > be a couple of large integral parameters. We divide 
/jy into K = 2Sm equal sub intervals, so that 



In J 



27rm 27rm 



for 1 < A; < ET. 

We then write the LHS of (ED as 



K 

(63) Var(Z+ - Z^^'+) = Cov{Znm - Z^M, Znm - Z^m) ■ 

ki,k2=l 

We divide the last summation into 3 different ranges. That is, we define 

(64) E^= ^ 

|fel-fc2|<l 
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(65) E2= ^ 

2<\ki-k2\<SR 

and 

K 

(66) ^3= ^ , 

\ki-k2\>SR 



and prove that for 1 < i < 3 



Ei 

lim — = 0. 

N^oo N 



5.2. Preliminaries. In addition to the covariance functions r = rj^ and 
j,M _ q£ and Y^^ respectively, defined on we introduce the joint 
covariance function 

(67) r*^'0(x) = 4^'°(x) = E [yiv(y)l^/(y + ^)] , 

which is a function of x indeed, by stationarity. Similarly to (jlSp . one has 
|^A^,0| < 1 by the Cauchy-Schwartz inequality. 

Using the spectral form (fSOt) (resp. ([56]) ) of Y/v (resp. Yp^^), one may 
compute the Fourier expansion of r^^'^ to be 



M,0(^) ^ v^(o):^^(o) + V 7^^;;^^ • ^ cos f ^x") 

n=l ^ ^ 

= a/ f(n)r^'^{n) ■ ^ cos T — x\ 
^-^ V J 'Km \m J 

n=l ^ ^ ^ 

1 ^ / ^ 



n=l ^ ' 



n 

cos I — X 

m 



In particular, rjJlf''^ is even, and 
(69) 

rN'^i'>T') = \/rN{n)r%^{n) = (^^^ J^-Xi<\n\<N{n)-^J (xi<|n|<JV * SM){n). 

Recall that to determine the second moment of a process X, we 

naturally encountered the distribution of the random vector ()16p . Similarly, 
to evaluate the covariances in ()63p , one naturally encounters the distribution 
of vectors 

= {Y;f{x,),Y^{x2),Yf{x^),Yf{x2)) 
with probability density (j)^'^ {ui,U2,vi,V2) and 

W2 = {YN{xi),Y^'{x2),Y;j{xi),Yf{x2)) 



with probability density (j)^'^Q{ui,U2,vi,V2) for some xi,X2 € In- As be- 

d (p^i^'^Q depend only on x = X2 

^N,M — 9n,M '^"Cl 9n,M,0 ~ ^N,M,0 



fore, the distributions 4h^'^ and (p^fJ'fln depend only on x = 0:2 — a^i by 



stationarity, and we will denote (p% = 4''n'm 4'% mo ~ 'Pl^'tfo accord- 
ingly. 
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Similarly to the mean zero Gaussian distribution with covariance matrix 
p!7|) of the random vector (fT6]) . both Wi and W2 are mean zero Gaussian 
with covariance matrices 



(70) Sw,m(x) 



/ 1 







rfj'ix) 


\ 




1 












-r^ {x) 


\M ' 


-rfl"{x 









-rf'ix) 


\M ' 


/ 



and 



(71) T^N^Mflix) 



1 



Yn Kx) 



Mflf X 



A/,0'/ N 
-r^' (x) 



-r^' (x) 
X' 

^2,Af 
-r^' (x) 



> ( 





-r^' (x) 



where, as usual, we denote 

A'2,;v := --^^(0); := -rr(0). 

Similarly to S(t) in (fT7|l . both Tjn,m{x) and SAr^M,o(2;) are nonsingular 
for X 7^ 0, and so 



(72) 

and 
(73) 



1 



(lir)'^ det T,N,Mia 



(27r)2y/detSAr,M,o(a;) 



=e 2 



We denote 
(74) 

and define the random vector 

(Fi = v^{x),V2 = V2(x)) = (y^(o),y/'(x)) 

conditioned upon Y]\f{0) = Y^^ (x) = with probability density function 
V'at M o(^i' ^2)- The random vector {Vi,V2) has a mean zero Gaussian dis- 
tribution with covariance matrix 
(75) 



( 



A: 



2-^ l-r-^-«(x)2 



Af,0" 



(x) 



'AT (2;) 
A2,Af 



N 



The matrix f^^^jvfo regular for x / 0. We have, analogously to 



'^Ar,Af,o(^^l>^^2) 



V'Af,A/,o(^l'^2) 
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Similarly, let jv/(^i' ''^s) be the probability density function of (Y^^'(O), Y^^' 



conditioned upon {0) = Y^{x) = 0. One then has 



(76) </'Af,M(^^l,^^2) 



27r^l-<(x)2' 



5.3. Auxiliary Lemmas. 

Lemma 5.1. One has the following estimates 
(1) 

O 



(2) 



(3) 



(4) 



(5) 



(6) 



I M II 
FtV -^iv||L2(%) 



I -'^.O II ^ 



MV4 



I M" // II ,01-'^ 



V -rivllL2(/^) 



MV4 



I A^' 'II ^ \ 

K -r^\\LHi.)=0[-^), 



M,0 

'•n 



-^'N\\LHi^)=o(^j^y 



Proof. First, we notice that (resp. ([U])) follows from ([T]) with ([2]) (resp. 
([2]) with dH)) by integration by parts and the Cauchy-Schwartz inequality. 

By the symmetry of all the functions involved, it is sufficient to bound 
II • ||^2(/+)- To establish ([1]), we note that for |x| < M, one has 



(77) Sm{x) = 1 + 



X \ 
M ' 



and both r^r and are rapidly decaying for large x, since 

« i 

X 
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and Sm is bounded, with constants independent of N or M. Thus, we have 

vrm 

Ikiv -^Af|li2(7^) = j {rfJ{x)-rNix)fdx 





M ■Km 



{rN{x){l - Sm{x))) 'i^'^J^ J "^Nix) X dx+ J r^ix) dx 

M 

M irm M 

« ^ jMxfx^) ■ x^dx + / 5 « ^ /-^'ix + ^ « ^. 



Af 



It is easy to estabhsh ^ using a similar approach. 

To prove ([2]), we wih use the Fourier series representation (fl9|) of r^v, and 
its analogue ([68]) for r^'^ with Parseval's identity. We then have by ([69|) 

II M,0 ||2 _ II -Af,0 - ||2 

Fat ~ ^A^IIl2(/jv) ~ ll'^Af ~''A^II/2(Z) 

N N 

= 2Y,rN{n) ■ WrN{n)- ^f^^{n)f < 2^f^(n) • jr^vW - fjlf (n)|, 

n=l n=l 

since for a, 6 > 0, 

(78) [a-bf <\a^ 

Continuing, we use the Cauchy- Schwartz inequality to obtain 

II A//,0 ||2 II - II II - 'Af II 

\Vn - '^NWi^i^l^) < lrAf||z2(Z) • IFA^ - ^A^ II«2(Z) 

Mr 1 



FAf||i2(2) • llrAT - 11,2(2) < 



M 



by (HI) of the present Lemma, and the obvious estimate ||r7v|| ^ 1- This 
proves part ([2]) of this Lemma. 

It is now easy to establish part @ of the present Lemma, using 

f"{n) = ^f{n). 

□ 

Lemma 5.2. The functions rjv(x), rjj{x), r^'^{x) and their first couple of 
derivatives are Lipschitz, uniformly with N, i.e. satisfy 

(79) \h{x) - h{y)\ < A\x - y\ 

for some universal constant A > 0. 

N 

Proof. The statement is clear for riv(x) = ^^^(m^)' ^^^^ ^Ni^) ~ 

n=l 

r]\f{x)SM{x) (due to the fact that Sm and S'j^j are bounded). 

It then remains to prove the result for r^^'^. From the representation (|68p 
it is clear that it would be sufficient to prove that the coefficients 



DISTRIBUTION OF ZEROS 



29 



are uniformly bounded. We will bound the square 

J^:;^-iXl<\n\<N*SM){n). 



(80) 



Using with ([55]) . we bound Siyiin) by 



(81) 



Suin) < < 



,15/2 / Sin 



M 



nM \ 
m } 



n = 



so that the coefficients ([80]) are bounded by 



m 



15/2 



n+7V 

E 

k=n-N 



sm 



fcAf A \ 8 



(MM) 



+ 



M 



m \/m 



M v-^ 



sm 



<kM\\^ N'T 



kM j ' M'^ 

m 

1 



N' I M 



N M M'^ {N/my 



+ K 1. 



This proves that the squared coefficients (j80|) are uniformly bounded, and 
thus, that rjl) 
stant ^4 > 0. 



thus, that r^'^ satisfy the Lipshitz condition (j79p with some universal con 



Lemma 5.3. Let I = [a, b] he any interval and h : I 
function satisfying ([75]) . Then for every x £ I, 

\h{x)\ < 2Ay^h\\%l^, 

provided that 



□ 

a Lipschitz 



(82) 

Proof. Let xq £ I and 

J := /n 



b - a> 



max /i(x) 



2y4 

|/i(xo)| 



2A 



\Kxo)\ 
2A 



Our assumption (j82p implies that interval J has length |J| > ^^^'^'^^ , and 
moreover, on J we have 



by (|79|) . Thus we have 



2 

L2(/) 



> / h{xYdx > \ J\ 



hixof ^ \hixo)\ h{xof 



> 



2A 



8A 



which is equivalent to the statement of this Lemma. 
Lemmas 15.11 15.21 and 15.31 together imply 



n 
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Corollary 5.4. For every x E Im, one has 

\rflix)-rN{x)\=0 



kiv'°(a^) -rN{x)\ = 



1 
1 

MV3 /' 



\rf{x)-r'Ax)\ = 0(^-^^ 
\r^'^'{x)-r'^{x)\=o(^j^ 
\rr\x)-rUx)\=0^-^ 
\rf'"ix)-r';,ix)\=o( ' 



. MV6 J ' 

uniformly w.r.t. x and N. 
Lemma 5.5. For every x £ I, we have 



E(y/(x)-yMx))' = o(-^ 



with the constant involved in the O -notation universal. 

Proof. By the stationarity we may assume that x = 0. We have by (jSOp and 

m 



E(y^(o) - y^(o))' = ^4^(0) + ^ f; (V^) - 

^ ^ n=l ^ ' 



Since r^ is supported in n < A^, we have 
(83) 



^ n=0 ^ /V „>2Ar 

We use ([78]) again to bound the first summation of ([83|) . getting 



2Af /t: 2Ar 



'■Km — ' x'TTm 

n=0 n=l 



, / 2Af . 1/2 

^ ^ n=l ^ 



by the Cauchy-Schwartz inequality, Parseval's identity and Lemma [5. H part 
©• 

To bound the second summation in (|83p . we reuse the estimate (j81|) to 
obtain 

iVl5/2 1 

SAiin) « • ^ ^ 
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SO that (|53|) implies 

N N 



k=-N k=~N ^ ' 

^13/2 ^ ^15/2 ^ 

and thus the second summation in ([83]) is bounded by 

^-'^ n>2N 

This concludes the proof of this lemma. 



□ 



Lemma 5.6 (Cuzick |CZj . lemma 4). Let Vi and V2 he a mean zero Gaussian 
pair of random variables and let 

n nr KV1V2-EV1EV2 
p = CoriV,,V2) -^7^p^=^- 

Then 

0<Cori\Vil\V2\)<p\ 

Remark 5.7. Lemma 15.61 mav be also obtained computing explicitly both 
sides of the inequality using the integral ()26p due to Bleher and Di [BD]. 

5.4. Proof of Proposition l473l 

Proof of Proposition Recall that (resp. Z^) is the number of the 
zeros of Yn (resp. Y^^) on In = [— 7rm,7rm]. The process Y^ is given 
in its spectral form ()56p with the RHS absolutely convergent, uniformly 
w.r.t. X € In- The rapid decay of implies that Y^^ is almost surely 
continuously differentiable, and we may differentiate (|56|) term by term. 

As stated before, for notational convenience, rather than showing the 
original statement of the Proposition, we are going to prove (j6ip . We want 
to bound Ei defined by ([M|) - ([66|) given the large integral parameters 5 and 
R. 

5.4.1. Bounding Ei. For every x G In, let Xn (I'esp. Xn m) indicator 
of the event {yiv(0)YAr(x) < 0} (resp. {Y^{0)Y^{x) < 0}). Intuitively, for 
S large (i.e. In,i is short) one expects at most one zero of either Yn or Y^ 
on In,i- Thus the number of zeros of Yn (resp. Y^^) on In,i = [0, r] with 
T := jg should, with high probability, equal x = X'n (I'esp. x*^ = X'n m)- 
Using the Cauchy- Schwartz inequality, we have for every ki,k2 

Cov{ZN,ki - Zn^u-^^Zn^ - ^ 



< ^Y^t{Znm-Z^'m) ■ \I^^<ZNM-Zf,\k,) = Var(Z^,i - Zj^i) 
by the stationarity. Therefore, 

El < SN\ai{ZN,i - Zf}^i) < SN¥.{Zn,i - Z^^^f 

^^^^ « sn( e{Zn,i - X? + nx - x"'? + nx'' - Z%'A . 
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We recognize the second summand of (j84p as the probabihty Pr{x 7^ X^^)- 
We may bound it as 

E(x - X*')' = Pr{x + X*') < Pr(sgn(y^(0)) / sgn(y/(0))) 
+ Pr(sgn(y^(r)) / sgn(F/(r))). 

We bound the first summand of the RHS of ()85p . and similarly the second 

one. For every e > 0, we have 

(86) 

Pr(sgn(y,v(0)) / sgn(y^^(0))) < Pr(|y,v(0)| < e)+Pr(| Yiv(0)-y/(0)| > e). 
The first summand of (j86p is bounded by 

Pr(|y;v(0)| <e) = 0(e), 
since YAf(O) is (0, 1)-Gaussian, and the second one is 

Pr(|y^(o)-y*''(o)|>6)«-^, 

by Lemma 15.51 and Chebyshev's inequality. 
Hence, we obtain the bound 

Pr(sgn(yjv(0)) / sgn(y^'^(0))) = ofe + ^ 



and, similarly, 

PT{mhoxsgn{Yn(T)) / sgn(y^^(r))) = o('e+ -A_Y 



Plugging the last couple of estimates into ([55]) yields that for every e > 

(87) EU-x")^=0(. + -^). 

The RHS of (j87p can be made arbitrarily small. 

Now we treat the third summand of (j84p . and similarly, but easier, the 
first one. We have 



k=l 

oo 

A:=2 

and and using the same approach as in (|57p in addition to some easy ma- 
nipulations yields 

r 

E[(Z^i)2 - Z^i] ^ j{r-x). k{x)dx, 


recalling the notation t := where 



^(x) = -:ryr (Vl -/32 + /)arcsinp) 

(1 — r^j ' 
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with notations as in (|58p . We saw already that K{x) is bounded, uniformly 
w.r.t. A^, so that 



Plugging the last estimate into (j88|) . we obtain the bound 

52 



(89) nz^,-x^r = o{i, 



and similarly, 

1 



(90) E(Zjv,i-x) =0. ^2 

as well. 

Collecting the bounds for various summands of (|84p we encountered i.e. 
(HZD, ([89]) and dlOD, we obtain the bound 



or, equivalently, 

which could be made arbitrarily small. 

5.4.2. Bounding E2. We write E2 as 
(92) 

^2 = '^Cov{ZN,ki - Z^^^^ZnM - ZN^ki) = '^^ZN,ki ■ {ZnM - ^A^fca) 
- y EZj^^fc^ • {ZN,k2 - ^A^fca) ~ [Znm - Zfl^ki] ^ [Znm - Z^^ki] 

and bound each of the summands of (|92p separately. In fact, we will only 
bound the contribution of the summands £'2,* and of the (slightly more 
difficult of the remaining two) -£2,21 bounding £^2,1 hi a similar manner. 
Note that the number of summands in each of the summations in (j92p . 
which equals the number of pairs {ki,k2) with 2 < |A;i — ^2! < RS, is of 
order K ■ SR = NS ■ SR = NS^R. 

We reuse the notation t '■= Note that differentiating ()52p yields 

rr(o)=-^^(o)+o(^, 

using ()77p near the origin. One then has 

^[Znm - Z^^ki\HZN,k2 - Z^\k.} = {E.[Zn,i - Zfji\f 

«{r{rUO)-rf'm)'«^„ 
by the stationarity, formula ([6]) and its analogue for Y^^ . Therefore 
(93) E2,t « N^. 
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Note that for |A;i — A;2| > 2, the intervals Inm and Inm are disjoint. 
Using a similar approach to |CLj , we find that 

(^2— fcl+l)T 



(A;2 — fci — 1)t 



\vi\ ■ \v2\{4>N,M,oi'"l:V2) - <^Ar,M (^^1 . ^^2 )) d^^l^f 2 



dx, 



using the notations ([7^) . We then bound £'2,2 as 
(94) 

E2,2 < NS^R . max {EZ*^,^ • (Z^,,, - 

|fei-fe2|>2 

■^NR max <^ \vi\\v2\ ■ \4>n m - 4>n Mivi,V2)\dvidv2} , 

T<X<TTm / / II, 



where we used the obvious inequality 

max < max 

2<|fci-fc2|<-R5 |fci-fc2|>2 

To bound the last integral, we exploit the fact that on any compact subset 
of M? we have 

l'^Af,A/,o(^^l'^^2) - (f>N,Mivi,V2)\ ^ 

as N ^ oo, uniformly w.r.t. x > t, whereas outside both (/)% M o ^^'^ 'Pnm 
are rapidly decaying. More precisely, let T > be a large parameter. We 
write 



(95) JJ = JJ + J J =:J, + J2. 

IR2 [-r,T]2 max{|t)i|}>T 

While bounding J2, we may assume with no loss of generality, that 

\vi\ > T 

on the domain of the integration. Let 



J2,l ■=JJ I'^llb2| • \(l)N^M,o{vi,V2)\dVidv2 
\vi\>T 

and 

^2,2 ■= jj l^lll^2| • \(f>N,Mi'^l^V2)\dvidV2, 
\vi\>T 

SO that 

(96) J2 < J2,l + J2,2- 
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Upon using the Cauchy-Schwartz inequality, we obtain 
(97) 



J2,2 < / dv2 J \vi\\v2\ ■ ^N^Mivi,V2)dvi 
\vi\>T 



1/2 / °? \ 1/2 



< ( / vldv2 / (f>% M{vi,V2)dvi] • dV2 vl4)% j^j{vi,V2)dvi 



T -oo T 

oo oo 



by dZel). 

Computing explicitly, we have 

(98) ^[Y^'\xf\Yii\^) = y/(x) = 0] = - '^''if = 0(1), 

where A^^' := —r^"{0), and, changing the order of integration, 

oo oo oo ^ 

-oo ^ i ^ ^vra^l-r*f(x)2 

where 



:= E[yi^'^'(o)2|y/(o) = y*^(x) = 0] = Afy - 



AT 

as well. Continuing, we bound the integral by 

oo oo 



X)2 



,2 



(99) 



-oo T 



dv2 / — = I v''^ exp ( — -v''^ ) dv' 



2 



« ^2 • - « 



l-rj^^(x)2 r2Jl-r]^(x)2 



(say), by the rapid decay of the exponential. 

Plugging (f98|) and (j99|) into (j97|) . and using the crude estimate 

l-r*f(x)»T2 

for T < X < vrm, we obtain the estimate 



(100) J2,2 « ^ , ^ • ^ « 
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Repeating all of the above for J2,i, we obtain 



(101) J2,l < 



1/2 



Mfi, 

'■n 



1 S 



using the same estimate 

1 - rf\x) » r^ 

which is easy to obtain using (|68p and Lemmas 15.11 and 15.21 

Plugging the inequality (jlOip together with (jlOOp into (|96|) . we obtain 



(102) 



J2 = 



Now we are going to bound Ji. Recall the definition ()74p of 
6^ Af 'with ([72]) and ([75]) , and the covariance matrices (f70|) , (ITT]) . 
Corollary 15.41 implies 



and 



o 



1 



MV6 



(here and anywhere else the inequality M < y where M is a matrix and 
2/ is a number means that all the entries of M are < than y). Expanding 
the determinants det Sjif and detS^''' into Taylor polynomial around the 
origin shows that they are bounded away from zero in the sense that for 
T < X < 7rm, 



detSjif, detS 



M,0 
N 



> r"" > 



1 

5^ 



for some constant A > 0. 



Remark 5.8. An explicit computation shows that det T,}^ (x) ^ x and 



also detT,'^'^{x) ^ x^, with universal constants. 



Thus also 



N ) 



M,Ox-l| 
N > I 



o 



and 



det S*^ 



detS 



N 



« o 



Ml/6 ) 

( 



MV6 



Substituting the estimates above into ([72]) and ([73]). and using ([71]) . we 
obtain 

g2A r^2 S'^^T'^ 

uniformly for t < x < Tim and \vi\ < T, where we used the trivial estimate 
|e^ — e^l < |rc — y| for x,y < 0. 

Integrating the last estimate for \vi\ <T and substituting into the defini- 
tion of Ji, we finally obtain 



(103) 



Ji = 



MV6 
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Upon combining ()103p and (jl02p . and recalling with ([95]) . we finally 
obtain a bound for E2 2 



(104) 



+ o|f 



and repeating the same computation for -E'2,1, we may find that the same 
bound is applicable for £^2,1 



(105) 



E2i = NR - [O 



Ml/6 



+ 



Using (|1U5|) together with (|1U4|) and (^5]) . and noting 
obtain a bound for E2 



we finally 



E2 = NR - [O 



V Ml/6 



+ 



fRS_ 



so that 
(106) 

which could be made arbitrarily small. 
5.4.3. Bounding E3. By the symmetry 

we may rewrite E^, as 
(107) 



M ^ 



E' 



3,1 



2-E3.2 + -E': 



3,3- 



First we treat the "mixed" term -£3,2) providing a similar treatment for 
the other terms. Assume with no loss of generality, that k2 > ki. Here we 
employ the random vector (Vi, V2) defined in section [5^21 Using the theory 
developed in [CLJ . modified to treat the covariance (see also remark [2^ . we 
may write 



(fc2-fcl+l)T 



Cov{Znm^Z^^^^) = ^ j 



(fc2 — fci — 1)t 

^[\Vi{x)V2{x)\\ 



{t -\x - {k2 - ki)T\) X 
-E|yj^(0)|E|y/'(x 



dx. 



where, as usual, we denote t := jg. Summing that up for \k2 — ki\ > SR, 
and using the stationarity, we obtain the bound 



-^3,2 < N 



E[\Viix)V2ix)\] 



E|y^(o)|E|y^'(x 



TVR 

2 



1 - r^'°ix) 



dx, 



38 



ANDREW GRANVILLE AND IGOR WIGMAN 



SO that 



(108) 



N 



2 



5Mlf§fa-E|r;(o,|E|rrwi 



dx. 



To bound the mtegral on the RHS of (|lU8p . we use the triangle inequahty 
to write 



IT 



|Cot;(|Vi|, 1^21)1 



(109) 



■kR 
2 



dx 



+ 



ttR 
2 



E|y^(o)|E|y/'(x 



1 - r^' {xf 



dx =: Js,! + J3,2. 



For ^ < x < nm, R sufficiently large, r^'^ is bounded away from 1 (see 
Corollary 15. 4p . and therefore, while bounding J^^i, we may disregard the 
denominator of the first integrand in (jlOOp . Note that if 1/ is a mean zero 
Gaussian random variable, then 



(110) 

and 

(111) 



Ei\V\) 



Var(|y|) 



yi-TVa^. 



1 - jVar(y). 



Note also that for ^ < a; < irm the variances Var(Vi(x)) and Var(V2(x)), 
given by the diagonal entries of ()75p . are bounded away from 0. This follows 



from the decay of r^'^{x) and r^'^ (x) for large values of x, due to Corollary 
15.41 Thus, an application of Lemma EU] yields 



< Cov{\V,l\V2\) < ^i=i^£^ffi^ « CoviVu V2f. 
^/YaT{Vi) ■ Var(V2) 



All in all, we obtain the estimate 



J3,l < J Cov{Vi,V2)'' 



dx, 



ttV 
2 
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which (this time, using the off-diagonal elements of (j75p ) is 
(112) 
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-^3,1 <■ 



M,0/ \ iviM' / \ 

Mfii't . r^' {x) ■ r^' (x) 



M,0', \2-\2 



N 



(x) 



2 



1 - rjlf'°(x)2 



dx <^ (^7v'° (^)^ + '^iv'^C^;) )dx 



2 



< W^n'^ - ''Jv|li2(^^) + llr^''''" - r'NWh^i^) + I {r'Nixf + rN{xf)dx 



M,0" // ||2 



2 



1 1 



by the triangle inequality, Lemma |5.H and the decay 



rN{x),r%{x) < -. 



To bound J3.2, we note that for ^ < x < irm, we may expand 



l + 0{r'J'\xf) 



Y^l -rj|f'°(2;)2 



with the constant involved in the 'O'-notation being uniform, since r^'^ is 
bounded away from 1 (by Corollary [531 say). Thus we may use the triangle 
inequality to write 



(113) 



E\Vi\-E\V2\-'E\Yl^{0)\E\Y^'{x 



71 7 n 

dx+ j r^'^{xf ■E\Vi\E\V2\dx 



TVR 

2 



ttR 
2 



< 



71 IfL 71 tn 

J E\Vi\{\E\V2\-E\Yf{x)\\)dx+ J E\Yf{x)\{\E\Vi\-E\YN'm\)dx 



ttR 
2 



nR 
2 



71 in 

+ j r^'\xf ■E\Vi\E\V2\dx. 



ttR 
2 



Now, (fTTOjl allows us to compute E\Vi\, E\V2\, IE|y^(0)| and E\Y^'^' {x)\; 
()75p implies that all of the expectations above are uniformly bounded for 
^ < a; < vrm. Thus the third term of ()113p is bounded by 



« / r'J^\xfdx « 1 + 



1 1 



ttR 
2 
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as before. We bound the first summand of (jll3p as 



< / \E\V2\-K\Y^^'{x)\\dx<^ 

2 2 



A^^' - A, 



dx 



Trm q/ 7rm 

« ^— f ''^^ « / r"-<>'lx?dx « — + i 



2 2 

as earlier, since A2^^' is bounded away from 0, and r^'^ is bounded away 
from 1 on the domain of the integration. The second summand of ()113p is 
bounded similarly, resulting in the same bound. Therefore 

(114) + i 

Recalling (fTU^ . the estimates (fTT^ and (flTID imply 

« J_ + 1 

by (|1U9|) . Bounding E^^i and £"3^3 in a similar (but easier) way, we get (see 

(nnzi)) 

5.4.4. Collecting all the estimates. Collecting the estimates (|9T]l . (|106p and 
(|115p . we see that 

Var(Z+ - ^Ei E2 E3 

N N N N 

^/ 5 1\ ^[RS^^T^ RS\ ^( 1 1 

= 0\eS + + 7; + O + — + O ^= + - 

which could be made arbitrarily small, upon making an appropriate choice 
for the parameters e, 5, R and T. Proposition 14.31 is now proved. 

□ 

Appendix A. The third moment of Z^^ on short intervals is 

BOUNDED 

Proposition A.l. Let L he a constant, K = NL and for 1 < k < K 

let Z^^f^ be the number of zeros of Y^' on [{k - 1)^, k^]. Then for L 
sufficiently large, all the third moments E(Zj^^^)^ are uniformly hounded by 
a constant, independent of N and k. 

Proof. By stationarity, we may assume that k = 1. For any r > let 
Z = Z^{t) be the number of zeros of y = Y^^ on [0, r]. Since L is 
arbitrarily large, we may reduce the statement of the present Proposition to 
bounding KZ^ (t)^: for t > sufficiently small. It will be convenient to use 
the shortcut r = r¥ . 



DISTRIBUTION OF ZEROS 41 

Using the formula for the high combinatorial moments of the number of 
crossings of stationary processes [CL) (see remarks 12.31 and I2.4p , we obtain 
the bound 

E[Z(Z- l)(Z-2)] < jj P{x,y)dxdy 

[0,r]2 

for the third combinatorial moment (the number of triples) of Z{t), where 
P = P^^ is given by 

.^^g. p. . ^ E[\Y'iO)Y'ix)Y'{y)\\Y{0) = Y{x) = Yjy) = O] 

(2vr)3/2y7(^ 

with f{x, y) = f^^{x, y) = detA with 

/ 1 r(x) r(y) \ 

A = r{x) 1 r{y — x) \ . 
\r{y) r{y - x) 1 / 

It is easy to compute / explicitly as 

(117) f{x, y) = 1 — r{x)'^ — r{y)'^ — r{y — x)'^ + 2r{x)r{y)r{y — x). 
Since 

EZ = -^/^?^ = 0{l), 
vr 

and we proved that the second moments EZ^ are uniformly bounded, while 
proving Proposition 14.41 (see condition [2l of Theorem . it is then sufficient 
to prove that the function P(x, y) is uniformly bounded near the origin. 
Denote the random vector 

(Vi,F2,^3) = (>^'(o),r'(x),y'(y)) 

conditioned upon 1^(0) = Y{x) = Y{y) = 0. The random vector {Vi, V2, V3) 
has a mean zero multivariate Gaussian distribution and we have 

(118) = 

by the definition ()116p . 

Applying the Cauchy-Schwarts inequality twice implies the bound 

E[|Fiy2V^3|] < (Eyi2)l/2(]Ey4)l/4(]Ey4)l/4^ 

Let Vi = Vi{x, y) be the variance of Vi for i = 1,2, 3. Computing explicitly, 
we have Vi = j^^TZi{x, y), where 

7^l(x, y) := A'2 det A - r'{xf{l - r{yf) - r'{yf{l - r{xf) 
— 2r'{x)r'{y){r{x)r{y) — r{y — x)), 

'R-2{x, y) := A2 det A — r'{x)'^{l — r{y — x)^) — r'{x — y)^(l — r(x)^) 

— 2r'(x)r'(x — y){r{x)r{y — x) — r{y)), 

and 

7^3(x, y) := A'2 det A - r'{y)\l - r{y - xf) - r'{y - xf{l - r{yf) 

— 2r'{y)r'{y — x){r{y)r{y — x) — r(x)), 
where, as usual, we denote A2 := — r"(0). 
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We then have 

The uniform boundedness of P{x, y) around the origin (x, y) = (0, 0) then 
fohows from applying Lemmas lA. 21 and IA.4I on (|119p . bearing in mind (jllSp . 

□ 

Lemma A. 2. Let f{x,y) be defined by ()117p with r = . Then 

f{x,y) > x^y^(y - xf 

uniformly w.r.t. N in some (fixed) neighbourhood of the origin. 

Proof. Recall that r^{x) = rj\[{x) ■ Sm{x), and we assume that the neigh- 
bourhood is sufficiently small so that Sm is given by a single polynomial 
([5T]) of degree 7 in We assume with no loss of generality, that x,y > 0, 



M 

and furthermore, that y > x. Let 

e^{x) :=r*f(x)-l, 

and 

0m{x) := r^ix) - 1. 

Let also 

smx 

Ooo{x) = 1 

x 

be the limiting function. We will omit the parameters m and M, whenever 
there is no ambiguity. 

We rewrite the definition of f{x,y) as 

f{x,y) = -{eixf + 9iyf + eiy-xf) 

+ 2{e{x)e{y) + e{y)e{y -x) + e{x)e{y - x)) + 2e{x)e{y)e{y - x). 

It is easy to Taylor expand 9 = 9^{x) as 

9{x) = af^^x' + 0{x^), 
where the constant in the 'O'-notation is universal, and 

1 1 

where 02 = — ^ is the corresponding Taylor coefficient of the limiting func- 
tion 9oo. We rewrite it as 

so that 
(120) 

9{x)9{y)9{y-x) = -^xV{y-xf (^l + o(^x^ + y^ + {y-xf + ^ + ^ 

Thus, it remains to estimate 
/2(x, y) := 2{9{x)9{y)+9{y)9{y-x)+9{x)9{y-x))-{9{xf+9{yf+9{y-xf). 



a¥m = a2{l + 0{- + 
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Let 

oo 

(121) 0oo(a;) = ^a„a;'^, 

n=2 

where it is easy to compute a„ to be 



(122) a. 



(-1)" 

^ even 



otherwise 
Similarly, we expand 9m and 9^ into Taylor series 

oo 

and 



CO 

ra=2 



3M/^^ _ „M n 

/ ^ "'n,m'^ ■ 

n=2 



We need the following estimates concerning the Taylor coefficients of 

M 



and ^„ 



Lemma A. 3. (1) We have the following estimates for the coefficients 
of 9m, 



1 

a2n,m = «2n ( 1 + O ( ^ k 



m m? 



for n <^m and 



for n'^ m. We have a2n+i,m = for every n. 
(2) We have the following estimates for the coefficients of 9^ , 



<m = a2n( l + 0( - + ) 1+0 



1 n^W „/ 1 



m m? J J \Af2(2n — 5)! 



^ ' «2n+l,m - S _ 1 'i n > 

['^U/^(2n-6)!i' "-"^ 

/or n <^ m, and 



(124) <-«^"™(4^)"-'^^^ 
/or m. 

We postpone the proof of Lemma IA.3I until after the end of the proof of 
Lemma (|A.2p . 
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We write 



oo 



+ x\y - xy + x\y - xy - x(*+^') - y^'+^^ - {y - xf+^'^), 

adding the summands corresponding to and (j, «). We introduce the 
polynomials 

Fij{x,y) := x'y^ + x^y' + y\y - xy + y\y - xy + x\y - xy 
^^"^^^ + x^{y - xy - - - (y - xf+^^ G Z[x, y], 

so that 

oo 

and 



f2,oo{x,y)= ^ aiajFij{x,y). 



i,j=2 
i,j even 

Note that ^2^2 = and for every even tuple (z, j) 7^ (2, 2), 

- xy\Fij{x,y), 
so that in this case, we may define 

Hid{x,y) = Jy2^(^y'^l^y e z[x,y]. 

It is easy to compute -^2,4 and H4 2 to be 

H2,A{x,y) = H4^2{x,y) = -6. 
We claim that for \y\ < r, if r is sufficiently small, 
(126) \Hi,j{x,y)\ «max{2|x|,2|2/|}2(^+^-3) 

To prove our claim, we write 
F2i,2j{x,y) 



{y - xy 

y — X y — X 

The sum of the coefficients of the homogeneous monomials of the last polyno- 
mial is <C 2^(*+-'^^). Since dividing by x'^y'^ does not increase the coefficients, 
the Hi J are bounded by 

< 22(*+^-i) •max{|x|,|y|}2(^+^'-=^) < max{2|x|, 2|y|}2(*+J-3), 

which is our claim (jl26p . 
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We then have 

/2,oo(a;,y) 



120204+ ^ a2ia2jHij{x,y), 



so that on any fixed neighbourhood of the origin 



+ 120204 



< ^ \a2ia2j\\Hi^j{x,y)\ 



x'^y'^{y — xY 



i+j>4 



by (fT22]l and (fT26]l . 

Thus to finish the proof of Lemma lA.2l it is sufficient to bound 

/2,oo(a;,y) - fiLix,y) 



x^y'^{y — x)^ 



We have 

f2,oo{x,y) - f2in{x,y) 



x'^y'^iy — x)2 



where 



We have 



i+j>6 



Fi,jix,y) 



x^y^iy — xy 



^ etien _j_ j-i odd _j_ 2 misled 



jj' odd 



mixed 



I, J even 



E 

i odd, J even 



|S--|= 5] |a2t„.a^,„- 02,02,||/f,,(x,y)| 

i+i>3 

1 i"^ +f. 1 



^ ( ~^ 



+ 



m m2 '(2i + l)!(2j + 1)! M2 • (2j - 5)!(2i - 5) 



• r 



2(i+i-3) 



{2i-5)\ (47rm)2i 



47rm 

e ^•0(l)^2(i+j-3) 



(47rm)2(^+i) 



nO(1) 2(i+i-3) ^ , J_ 



which is sufficient. 

The main problem while treating the odd and the mixed terms is that 
for such a tuple Fij{x,y) is not necessarily divisible by x'^y'^{y — x)^. 

However, in any case, it is divisible by x'^{y — x)^. We recall then our 
assumption y > x, and thus 



x^y^{y — xY 



y2 x'^{y — x)2 



We then define for an odd i > 7 and any j (i.e. j > 2 and j 7^ 3, 5) the 
polynomial 



Gi,j{x,y) 



F 



x'^{y — x) 



2 ■ 
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The polynomial Gij is a homogeneous polynomial of degree i + j — 4 and 
we have 



(127) <max{2|x|,2|y|} 
similarly to ()126p . Thus, one has 
1 



i+j-i 



(128) 



< max{2|x|,2|y|} 



i+j- 



(here we use < |y|). 
We then write 



yodd _ \ " M M Gi^j{x, y) I 

using the same approach as in case of E'^*'^", this time plugging (|128p . Sim- 
ilarly, one obtains the estimate 



ymixed " , 

Combining ()120p with the estimates on shows that 
f[x,y) , r>f 2 , 2,1, 1 



1 1 



x^y^(y — x)^ \ "T- M^y 

where 

a = -120204 - 2/6^ = — . 

135 

This concludes the proof of Lemma IA.2I assuming Lemma IA.3I 



□ 



Proof of Lemma \A.3[ First, it is clear that part ([2]) of Lemma I A. 3 1 follows 
from part ([1]) by ()5ip . which holds for x > sufficiently small. 
We write 

''"■(^'-^(^'•(^)-^ 

where h{x) = ^jf-j- The multiplication by poses no problem here. 

Now, the Taylor expansion of h{x) is well-known to be 

, „2n 

Mx) = ^(-l)"+H2^"-2)i?2„^, 

n>0 ^ 

where Bn are the Bernoulli numbers. Recalling the Taylor expansion 

sinx ■r;^ • x^-^ 



we obtain, after a little rearranging. 



V X \2m) ) ^\ ^ \ 2n )^ ' {Imf'' ] {2k + 



Now, 



k>l ^ j,n>0 
j+n=k 



i?2.~(-l)"-^-2 ^^"^ 



(2vr) 



2?i ■ 



and therefore the coefficient of (2A:+i)! ^® 
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1 + E 



2A; + 1\ (2n)! 



\2n 



^ \ 2n ) (47rm) 

n>l ^ ^ ^ ^ 

By comparing consecutive terms in this series, we find that this is 1 + 
Oil? jrr?^ provided that k <^ m, and it is 



4e7rm 

if /c ^ m, using Stirhng's formula, since the maximal term occurs when 

2n w 2A; + 1 - Anm. 



Hence we have 
for n <^ m, and 

for m. 
Lemma A. 4. Let 



1 

a2n,m = «2n ( 1 + O ( ^ n 



1 \ 



47rm . 

□ 



7^l(x, y) := Xyix, y) - r'{xf{l - r{yf) - r'{yf{l - r{xf) 
— 2r' {x)r' {y){r{x)r{y) — r{y — x)), 

'R'2{x, y) := A2/(2;, y) - r'{xf{l - r{y - xf) - r'{x - yf{l - r{xf) 

— 2r'{x)r'{x — y){r{x)r{y — x) — r{y)), 

and 

'J^3ix, y) := A2/(x, y) - r'{yf{l - r{y - xf) - r'{y - xf{l - r{yf) 

— 2r'{y)r'{y — x){r{y)r{y — x) — r(x)), 

where r = rf^ and f is defined by ()117p . Then 

ni{x,y) = o{xY{y - xf), 

TZ2{x,y) = 0{xV{y-xf), 

and 

ns{x,y) = o{xY{y - xf). 

Proof. We prove the first statement only, the other ones being symmetrical. 
We will assume that x,y > and moreover y > x, so that Sm is given 
by a polynomial (j5ip of degree 7 in For brevity, we denote h{x,y) = 
h^{x,y) = TZi{x,y). Similarly, we denote hoo{x,y), defined the same way 
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as {x, y), where we use r(x) = roo(x) = rather than . We rewrite 
h{x, y) in terms of := r — 1 as 

h{x, y) := Xy{x, y) + 2 {e'{xfe{y) + 9'{yfe{x) 
-e'{x)9'{y){e{x) + e{y)-e{y-x))^ 

+ {e{x)e'{y)-e{y)e\x))\ 

Expanding 9 and 9' into the Taylor series as ear her and using = — 2a2, 
we have 
(129) 

hooix,y)= ^ a2jia2iia2i2a2j2^2ii,2ji,2i2,2i2(2;>2/) 

+ ^ a2ia2ja2kB2i,2j,2k{x,y) + ^ a2a2ia2ja2kC2i,2j,2k{x,y) 

i,j,k>2 jj>2,fe>l 

+ X] "'2"-2ia2jD2i,2jix,y) + ^ a2a2ia2jE2i,2jix,y) + '^ala2il2i{x, 

i,j>2 i,j>2 i>2 

+ '^ala2iJ2i{x,y), 

i>2 

where 

A,j,,i,j,ix,y) = hj2 ■ (x'^y^'-' - x^'-'y'^) ■ (x'^y^^-' - x^'-V), 
Bij,kix,y) = 2zj(x*+^--2/ + y^+J-^x>^ - x'-'y^-\x' + / -(y- , 

Cij,k = -ix'y^{y-x)^, 
Di,j{x,y) = -4(y - xy{x^y' + x'y^) + 4{yx' - xy'){yx^ - xy^) 
+ Ajix'y - xy'){x^y^-^ - x^'^y^), 

Eij{x, y) = -2F^j + Ai{2{x'y^ + x^y') - {x'-'y + f-'x){x^ + y^ - (y - x^)) 
+ 2ij(x*+^'-2y2 + 2/»+i-2^2 _ 22.iyj) ^ 

(here the polynomial Fij is defined as in p25|) ). 

Ii{x, y) = -4Fi,2 + 8(x2y^ + x^y^ _ xy{x^ + y' - {y - xf)) = 0, 
so that we may disregard li altogether, 

Ji{x, y) = -4x2y^(y - x)' - 4(y - xf{x^y' + x'y"^) 

+ 8xy{yx^ — xy*)(x — y) + 4i(x — y)(x'^y* — y^x*), 

and 

5 = z^(({(n,i)}x{(z2,i)})u({(n,i)}x{(i,i2)})u({(i,ii)}x{(f2,i)})). 

From all the above, it is easy to check that for all the (even) indexes within 
the frame, A,C,D and J are divisible by P{x,y) := x'^y'^{y — x)^, and, 
moreover, i?jj^fc(x,y) + Bj^i^k{x,y) and Eij{x,y) + Ej^i{x,y) are divisible 
by P{x,y) (in particular, -Bj^j^fc and -Bj^j are). It follows then that all the 
polynomials above vanish, unless their degree is > 10. In addition, we 
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have the following estimates, which follow from the same reasoning as while 
proving (jl27p : 



A, 



«1 Jl,«2 J2 



Bi,j,k{x,y) 



x*y^{y — x)^ 



'i;j,k 



Di,j{x,y) 



x^y'^iy — x)^ 



x'^y^iy — x)"^ 



x^y'^iy — x^ 



<max{2i+^|x|,2^+^|y|f , 



(say), for some e > 0, where d is the degree of the corresponding polynomial 
on the LHS. 
Thus, we have 



(130) 



hooix,y) 



x^y^iy — x)'^ 



0{x' + y^), 



by the rapid decay (jl22p of the Taylor coefficients of ^oo- The constant c 
may be computed explicitly to be 

1 



212625 

To conclude the proof of Lemma lA.41 in this case, we need to bound 

hoo{x,y) - h^{x,y) 



x'^y^iy — xY 
Similarly to ()129p . we have 

uM( x_ Sr^ I M „M M M A. . . . (j. 

«1 j'l,«2,j2GS' 

J- \" 'n^^ R I't- .A -L \" ' H /'o 

T" 2-^ "'i,m"'j,m"'k,m-^i,j,k{X, (J) ^ "'2,m^i,m^j,m^k,m'^i,j,k\-' 



«J>4,fc>2 



+ ('^2,m) ^i,m'^j,m-^ijix,y) + a2„ia^,^aj„^Eij{x,y) 

iJ>A i,j>4 

+ J2'Kmfarj^Mx,y), 



i>4 



where the ' in the stands for the fact that all of the indexes involved in 
the summations above are ^ 3, 5, and 



\ {(^1,2)} X {(i2,2)}U{(ii,2)} x{(2,j2)}U{(2,ii) x {(^2,2)} 



S 



As in the proof of Lemma [A. 21 the tricky part here is that for odd indexes, 
the polynomials are no longer divisible by P{x, y) = x^y'^{y — x)^. However, 
we notice that, by our assumption y > x, it is sufficient that they are 
divisible by x^(y — x)^, and their degree is > 10 (see the end of the proof 
of Lemma lA.2p . To check this, we note that there is certainly no problem 
with Ai^j^^i^j^{x,y), {Bij^k + Bj^i^k)ix,y) and Cij^k{x,y). Moreover, one 
can easily check that Dij and Eij + Ej^i are divisible by x^(y — x)^ for any 
z, j > 2, and for any i > 5, Ji is divisible by x^y'^{y — x)^. Finally, we check 
that all the relevant polynomials have degree > 10. 
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This, together with the rapid decay (jl23p for n <^ m and (|124p for n ^ m, 
of the Taylor coefficients, imply that 



hoo{x,y) - h^{x,y) 



1 1 

^ — + T79- 



Combining the last estimate with p30p concludes the proof of the present 
Lemma. 

□ 
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